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Summary. — We review recent theoretical results for hole spins influenced by 
spin-orbit coupling and Coulomb interaction in two-dimensional quantum wells as 
well as the decoherence of single hole spins in quantum dots due to hyperfine in- 
teraction with surrounding nuclear spins. After reviewing the difl'erent forms of 
spin-orbit coupling that are relevant for electrons and heavy holes in III-V semi- 
conductor quantum wells, we illustrate the combined effect of spin-orbit coupling 
and Coulomb interactions for hole systems on spin-dependent quasiparticle group 
velocities. We further analyze spin-echo decay for a single hole spin in a nuclear-spin 
bath, demonstrating that this decoherence source can be controlled in these systems 
by entering a motional-averaging regime. Throughout this review, we emphasize 
physical effects that are unique to hole spins (rather than electrons) in nanoscale 
systems. 



Introduction 



Spin-related phenomena in semiconduetor nanostruetures have become increasingly 
important in many branches of condensed matter physics, including the study of spin 
currents [1, 2] and spin-Hall effects [3, 4], the search for Majorana fermions in solid- 
state systems [5, 6, 7, 8, 9, 10, 11], topological insulators [12, 13], and applications of 
quantum information/computation [14, 15, 16, 17, 18, 19]. Improved understanding of 
these phenomena could lead to new technologies based on our ability to control electron 
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spins. The ultimate demonstration of this control would be the precise manipulation of 
the quantum states of independent spins as well as interactions between them. Such a 
level of control for spins in semiconductor quantum dots may allow for a scalable im- 
plementation of quantum computing [14, 17, 18, 16, 19]. Additionally, semiconductor 
'spintronic' devices relying on spin-polarization and spin currents (rather than electric 
charges and charge transport) could lead to low-power high-density devices that would 
be unthinkable with conventional electronics [15, 1, 2]. A central proposal in spintronics 
is the Datta-Das 'spin transistor' [20]. In the Datta-Das transistor, a spin-polarized cur- 
rent would be modulated through the Rashba spin-orbit interaction of electron carriers. 
Although effects very similar to those proposed in the original paper by Datta and Das 
have finally been demonstrated [21], there is still much to be understood about, e.g., the 
role of Coulomb interactions in these devices (the subject of Sec. 3). Intense research 
efforts on quantum dots in III-V semiconductors have also led to remarkable progress in 
single-spin initialization, manipulation, and readout. All the basic ingredients for quan- 
tum computation (the so-called "DiVincenzo criteria" [22]) have been demonstrated in 
isolated (separate) experiments for spins in quantum dots. However, there remain many 
questions about the viability of combining these "nuts-and-bolts" in the presence of a 
random solid-state environment. These questions will be partially addressed in Sec. 4. 

In addition to nanoscale devices based on electron spins, there is a growing interest in 
the spin properties of valence-band holes in semiconductors. We will focus our attention 
in this Chapter on holes in III-V semiconductors (GaAs, InAs, InSb, etc.), which share 
a common band structure. Much of the interest in hole-spin systems originates from 
the inherently strong spin-orbit coupling of the valence band. This strong spin-orbit 
coupling introduces unique features for hole carriers arising from the orbital moment 
of the p-like valence band, relative to the s-like conduction band for electrons in these 
materials. This review will therefore illustrate the specific advantages and disadvantages 
of hole-spin systems over their electron counterparts. We will analyze these differences 
as they apply to both spin transport (Sec. 3) as well as single-spin manipulation and 
decoherence (Sec. 4). 

In Sec. 2 we give a brief review of the dominant spin-orbit terms for electrons and 
holes in two dimensions. Since local electrical gating is typically much easier to achieve 
experimentally than local magnetic fields, spin-orbit interactions represent an essential 
tool to address the spin degree of freedom at the nanoscale and a vast body of literature 
exists on the subject [23, 2]. Most studies related to spin-orbit coupling are based on 
single-particle models, neglecting (or ignoring) the effects of Coulomb interaction. On 
the other hand, the role of electron-electron interactions in electron liquids is a well- 
developed and active field of research in itself [24], although most studies in this field 
neglect or ignore effects of spin-orbit coupling. A natural question that arises is then 
how the properties of such two-dimensional systems with strong spin-orbit interactions 
arc modified by the presence of Coulomb repulsion, the subject of Sec. 3. The need for 
a better understanding of this problem is highlighted by recent experiments on dilute 
hole liquids, in which strong correlations and spin-orbit coupling effects coexist. In this 
regime, these hole liquids show a significant deviation from the conventional behavior of 
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two-dimensional electron liquids [25, 26]. 

Interaction effects can be computed accurately in the high-density limit by controlled 
perturbative methods (and some non-perturbative arguments have been put forward as 
well [27, 28]). In contrast, a quantitative understanding of the more correlated low- 
density regime including spin-orbit coupling is not yet possible. Various aspects of 
this problem have been examined in the theoretical literature, including quasiparticle 
properties (effective mass and lifetime) [29, 30, 31, 32, 33, 34, 28], screening and plas- 
mon excitations [35, 36, 37, 34], and more exotic collective excitations, including chiral 
spin waves [38]. Other works have considered ground-state properties including the 
repopulation of the spin bands [39] , the Hartree-Fock theory of broken-symmetry spin- 
polarized states [40, 41, 42], and low-density inhomogeneous phases [43]. Studies of the 
ground-state energy including correlation effects can be found in Rcfs. [44, 33, 45], while 
effects of interactions on the spin susceptibility with Rashba spin-orbit coupling have 
been studied in Refs. [46, 41, 34, 47, 48], where several types of non-analytic correc- 
tion have also been computed [47, 48]. We would also like to mention the recent surge 
of interest on topological insulators and schemes supporting localized Majorana modes 
[12, 13, 5, 6, 7, 8, 9, 10, 11]. While these systems are characterized by strong spin-orbit 
couplings, a proper understanding of interaction effects is an emerging and interesting 
problem in these fields as well (see, e.g., [49, 50], for Majorana modes). 

Part of the above mentioned work (that concerned with hole physics) is reviewed in 
Sec. 3 of this Chapter. In Sec. 3, we will pay particular attention to the behavior of the 
exchange-correlation energy and quasiparticle excitations by contrasting expected effects 
in electron and hole liquids. As we will discuss, the group velocities of quasiparticles at 
the Fermi energy have quite distinct features for hole systems in the presence of spin- 
orbit coupling. A particularly vivid demonstration of this difference is provided by a hole 
wavepacket moving ballistically in a two-dimensional system. As we discuss in Sec. 3, an 
initially spin-unpolarized hole wavepacket at the Fermi energy will spatially separate into 
its two spin components [28]. In contrast, electrons moving in the presence of conventional 
Rashba and Dresselhaus spin-orbit interactions will remain spin-unpolarized. 

While spin transport and manipulation are important for potential spintronic devices, 
quantum coherence and control are necessary for potential implementations of spin-based 
quantum computing, as well as more general quantum information processing tasks. Hole 
spins confined to semiconductor quantum dots have several potential advantages from 
the viewpoint of quantum coherence as well. The main difference between hole spins 
and electron spins can be traced once again to strong spin-orbit coupling. The result 
of this coupling is to lock the orbital angular momentum / = 1 to the spin angular 
momentum, s — 1/2. The states at the top of the valence band (the heavy holes) then 
transform like states with angular momentum ~ ±3/2 [51]. A primary decoherence 
source for electron spins is the Fermi contact hyperfine interaction between an electron 
spin and surrounding nuclear spins [52, 53, 54, 55, 56]. For holes, the difference in the 
microscopic spin-orbital state drastically modifies its coupling to the nuclear spin bath, 
and consequently, its decoherence properties [57, 58]. Heavy- hole spins in quantum dots 
are only affected by a weaker and strongly anisotropic hyperfine interaction [57], leading 
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to longer coherence times in a transverse magnetic field [57] or for a "narrowed" nuclear- 
spin bath state [52, 53]. 

Although single- electron spin coherence measurements arc becoming more common, 
single- /loZe spin coherence is a relatively new frontier with many as- yet unanswered ques- 
tions. Recent experiments have determined the value [59, 60] and sign [61] of the hyper- 
fine coupling constants for hole spins and have achieved coherent optical control of the 
hole spin in single [62, 63] and double [64] quantum dots, thus enabling measurements 
of the relaxation time Ti [65, 66, 67], estimates of a bound on the ensemble-averaged 
free-induction decay time T2* [68] , and more recently, real-time measurements of Tj* 
the spin-echo decay time, T2 [62]. In Sec. 4 we will focus on a review of recent theoretical 
work exploring related hole-spin-echo decay in experimentally accessible regimes [69] . 

2. Types of spin-orbit coupling 

Spin-orbit couplings in semiconductors are discussed in several excellent books [23]. 
Here we will restrict ourselves to a brief outline of the main spin-orbit mechanisms at 
play in III-V quantum wells, by emphasizing an important difference between electron 
and hole systems: while for electrons the dominant spin-orbit coupling mechanisms are 
linear in momentum, the dominant spin-orbit terms for heavy holes have a non-linear 
(quadratic or cubic) dependence on momentum. This dependence has been confirmed 
in the literature on heavy holes, both theoretically and experimentally. As we will see 
in Section 3, this distinction is important in discussing the interplay of spin-orbit and 
Coulomb interactions in a two-dimensional liquid. Quantities that are modified by this 
interplay include the exchange-correlation energy and the effective mass and lifetime of 
quasiparticle excitations. 

2'1. Rashba and Dresselhaus spin-orbit interactions. - For two-dimensional electrons 
(with motion confined to the x — y plane) the two main contributions to the spin-orbit 
interaction arc well known. Adding these to the kinetic energy, 

(1) ^ ^ 9^ ^ OliVyC'x - PxO-y) + P{Py(7y - PxCTx), 

where p is the momentum operator, m the band mass, and cr is the vector of Pauli 
matrices. The first spin-orbit term (proportional to a) is the Rashba spin-orbit cou- 
pling [70, 71]. This term has the same structure as the Dirac spin-orbit coupling 
oc Vy(r) • (cr X p), if the potential ^(r) describes an electric field EzZ. In fact, this 
type of spin-orbit coupling is caused by an asymmetry in the confinement potential 
(along z) which defines the quantum well. To derive this contribution, one can apply 
perturbation theory to a multiband Hamiltonian for the envelope functions. This is the 
usual procedure applied to the 8-band Kane model (see, for example [23, 72]). This 
procedure yields an estimate of the coupling a as follows 
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where we use the notation of Rcf. [23] for the band structure parameters: Eq is the main 
gap of the III-V semiconductor, Ao is the spin-orbit sphtting at the top of the valence 
band, and moP/ft is a matrix element of the momentum operator between the valence 
and conduction bands (with mg the free-electron mass). Note that a would be zero for 
Aq = 0. For III-V semiconductors, P ^ 1 eV-nm, which gives Ha ~ 0.05 meV-nm for 
GaAs and Ha ^ 15 mcV-nm for InAs, using a value of the electric field ~1 V/jim. 

The second spin-orbit term in Eq. (1) (that proportional to 13) is the Dresselhaus spin- 
orbit coupling [73, 74]. This term arises from the bulk asymmetry of III-V semiconductors 
and for this reason it is already present in the three-dimensional band structure with the 
following form [73]: 

(3) Hd = 1d[{pI~pI)PxOx + {pI -pl)Py(Ty + {pI ~pl)p,a^], 

where x, y, z are along the crystallographic axes. An estimate of can be obtained by 
applying perturbation theory to an extended Kane model (including 14 bands), which 
leads to a formula analogous to Eq. (2). This procedure yields U^^d — 27 meV-nm^ 
for both GaAs and InAs [23]. Although Eq. (3) has been known for a very long time, 
there is some recent controversy related to the accuracy of this value for in GaAs. 
A significantly smaller value h?^D — 9 meV-nm^ was obtained in Ref. [75], from a 
theoretical analysis of conductance measurements through open quantum dots [76, 77]. 
A collection of theoretical and experimental values for H^^d can also be found in Ref. [75] . 
On the other hand, good agreement with spin-relaxation measurements in the impurity 
band of GaAs has been obtained in Ref. [78], using h?"jD — 27 meV-nm^. 

Starting from Eq. (3), we can see that the two-dimensional Dresselhaus coupling in 
Eq. (1) can be derived by substituting by the expectation values (p^) = and (p^) ~ 
{ttK/WY , where W is the width of the quantum well, and neglecting small terms that arc 
cubic in the in-plane momentum. Using h'^^D = 27 meV-nm"^ and = 20 nm gives hj3 ~ 
S^7_D(7r/VF)^ — 0.7 meV-nm^. So, with these specific values of Ez and W, we find that 
for GaAs the Dresselhaus spin-orbit coupling is dominant over the Rashba mechanism, 
while for InAs the opposite is true. This derivation also emphasizes that Eq. (1) only 
applies if the growth direction of the quantum well is [001], i.e., along a crystallographic 
direction [74]. From Eq. (3) we can also derive the analogue of Eq. (1) for different 
growth direction of the quantum well. An interesting case (see, e.g., Ref. [79]) is the 
growth direction [110], for which the Dresselhaus coupling becomes simply proportional 
to (Jx + f^yi thus the spin is conserved in the direction perpendicular to the quantum well 
growth direction. 

Several new phenomena arise from the simultaneous presence of both Rashba and 
Dresselhaus couplings. These include an anisotropic conductivity [80], the formation of 
'electron-beams' [81], and a beating of Friedel oscillations [82], to mention just a few 
examples. In general, the presence of both Rashba and Dresselhaus terms results in a 
Hamiltonian which is anisotropic in the x — y plane. The special (symmetric) case a = (3 
results in a new conservation law [83]. One consequence of this conservation law is the 
existence of long-lived spin textures (the so-called 'spin helix') [84, 85, 86]. Determining 
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the precise values of a and (3 in any given sample is often challenging. The ratio a//3 can 
be measured through the spin-galvanic effect [87, 88], while it is also possible to extract 
the individual values of a and /3 independently using Faraday rotation to monitor the 
precession of electron spins as they move through a sample [89] . 

A large body of literature considers the effect of pure Rashba spin-orbit interaction 
(/3 = 0), which is isotropic in the x — y plane. To justify this assumption, one can notice 
that the pure Rashba Hamiltonian is equivalent to the case of pure Dresselhaus (a = 0), 
after performing a suitable rotation in spin space Gx^y — > 'yy,x- Thus, it is sufficient that 
one of the two spin-orbit mechanisms is dominant over the other. This assumption is 
often reasonable since, as we have seen, the relative strength of a and /3 depends on 
the specific choice of materials and the electronic confinement along z. For example, 
the Rashba spin-orbit coupling is absent in symmetric quantum wells. Furthermore, 
other material systems exist where the Rashba spin-orbit coupling is important and the 
Dresselhaus term might vanish. The presence of a linear-in-momentum spin splitting has 
been revealed in angle-resolved photoemission spectra from the Shockley surface state on 
Au(lll) [90]. In that case, spin-resolved photoemission experiments [91] directly reveal 
the spin structure of the electronic states, which are found to agree well with the Rashba 
eigenstates [see Eq. (5), below]. Furthermore, in the quasi- two-dimensional electron 
gas at the LaAlOs/SrTiOg interface [92], a tunable Rashba spin-orbit interaction can 
coexist with superconductivity [93, 94]. For this system, the existence of persistent spin 
oscillations was predicted in Ref. [95] . The presence of Rashba spin-orbit coupling also 
plays an important role in thin asymmetric ferromagnetic layers in contact with heavy 
elements, including Ft/Co/ AlO^ structures [96, 97, 98, 99]. 

2'2. Generalized spin-orbit coupling. - We can generalize the discussion of the last 
section by considering the following model Hamiltonian [39, 28]: 

"""^2^ + '^ 2 ' 

where p± ~ px i: ipy, a± = (Tx ± ioy. Introducing Eq. (4) is useful since for n = 1,2,3 
it represents different realizable spin-orbit coupling terms. In particular, the n = 1 case 
of Eq. (4) is equivalent to the Rashba term [oc a in Eq. (1)]. Alternatively, the n = 2, 3 
Hamiltonians arc appropriate, for example, for heavy holes in III-V semiconductor quan- 
tum wells. With n = 3, Eq. (4) has been shown to describe the main contribution to 
the spin-orbit coupling for holes in the presence of an asymmetric confinement poten- 
tial [100], in good agreement with magnetoresistance measurements [101]. The dominant 
^ p^ dependence can be related to the total angular momentum of the heavy holes. 
While electrons in the s-like conduction band have angular momentum ±1/2 (from the 
spin), the t / i states of the effective heavy- hole Hamiltonian (4) with n = 3 correspond 
to Jz = ±3/2. This is a consequence of the p-likc character of the valence band, since 
the orbital angular momentum combines with the spin into J ~ 3/2 states and the con- 
finement along z further splits the degeneracy between heavy holes (Jz = ±3/2) and 
light holes ( Jz = ±1/2). Recently, an anomalous behavior of polarization in magnetic 
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Fig. 1. - Chiral energy bands of a spin-orbit coupled system and parameters characterizing the 
strength of the spin-orbit coupling, g, the Coulomb iteraction, r^, and the chirality, x- These 
parameters are written in terms of the densities, n±, and Fermi wavevectors k± of the two chiral 
bands. 



focusing experiments for holes [102, 103] has been successfully interpreted on the basis 
of the n = 3 Hamiltonian (4). Finally, the coupling for n = 2 is present for heavy-hole 
systems with a finite in-plane magnetic field B and j oc B [23, 104, 39], consistent with 
considerations on time- reversal symmetry (i.e., the Hamiltonian must be invariant under 
the simultaneous change p — >■ — p, cr — )• — cr, B ~B). 

Diagonalizing Hq yields the following eigenstatcs, with the corresponding energies: 

where k is a wavevector in the x — y plane, 9k is the angle k makes with the x-axis, L is 
the linear size of the system, and ± labels the two chiral spin branches. We will generally 
assume 7 > 0, such that the " — " subband has lower energy. 

In the absence of spin splitting, the Fermi wavevector is kp = ^/STmJ, with the 
total density. However, for 7^0, the densities n± of the two non-degenerate spin 
subbands are different (with Us — + nJ). In this general case, it is useful to write the 
two distinct Fermi wavevectors k± (sec Fig. 1) as follows: 

(6) fc± ^ kp^/WrxV 

When < X < Ij there are two Fermi circles with radii k± and we have that x = (n_ — 
nj^)/ns. Thus, x coincides in this case with the chirality of the system. The parameter 
X is analogous to the fractional polarization of spin-polarized systems (although in this 
case the spin quantization axis depends on the direction of k and so there is no definite 
spin polarization when all states of the same |k| are filled uniformly for all directions in 
/c-space). 

To describe the non-interacting properties, and more specifically the value of x, we 
note that the Hamiltonian (4) can be characterized for all n = 1, 2, 3 by the ratio of the 
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spin-orbit energy splitting Eso = '2jh"'kp to the Fermi energy Ep = h^kp/2m. We thus 
define the dimensionless parameter: 

which aUows us to write the ground-state occupancy as 



n = l, 

(8) Xoig) = { 9, " = 2, 

In all three cases, the leading dependence at small g is simply given by xo — 3- 

In writing Eq. (8), we have assumed the most common situation of a rather small g, 
such that Xoig) < 1- The limit of large g requires a separate discussion for n = 1, due to 
the fact that E^{k) of Eq. (5) has its minimum at a finite wavcvcctor, k = mj/h, rather 
than at fc = 0. Thus, for n = 1 and sufficiently low density, the occupation is in the 
" — " subband and forms a ring in momentum space. In this case, Eq. (6) maintains its 
validity but x ^ 1 ^-i^d k± now determine the inner and outer radii (/c_ < mj/h < /e-|_) 
for the ring of occupied states. The precise condition for this regime is g > \/2, for which 
we can write Xoid) = 3^/4 + 1/.9^ > 1. Note that, since the chirality for a single occupied 
band is simply (n_ — n+)/ns = 1, the formula x = [n-- ~ n-^.)/ns in Fig. 1 is no longer 
valid when x > 1- It is then more appropriate to refer to x as a 'generalized chirality'. 
On the other hand, both spin subbands Ej_ {k) of Eq. (5) have minimum energy at fc = 
if the spin-orbit splitting is quadratic or cubic in momentum (for n = 2,3, respectively). 
For n — 2,3, and in the limit of small densities, the occupation of the two chiral bands 
is therefore always determined by two disks with radii given by k± in Eq. (6). 

3. Spin-orbit coupling and Coulomb interactions 

We now turn to many-body effects and discuss a few significant consequences of 
the interplay between spin-orbit coupling and Coulomb interaction. For the same non- 
interacting problem discussed in Sec. 2'2, we consider here the following Hamiltonian 
("=1,2,3), 



(9) g = y [ii+z/^-"-'+'^^+"--]+ly 



2 ^ erh 



where is the dielectric constant and standard terms associated with the neutralizing 
background have been omitted for simplicity [24] . The presence of the Coulomb repulsion 
introduces another dimensionless parameter in the problem in addition to g: 

(10) 
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where ub = erh^/me^ is the effective Bohr radius. The Wigner-Seitz radius rg charac- 
terizes the strength of the Coulomb interaction with respect to Ep: in the small-rs hmit 
(approached at large density Ug), the Coulomb interaction can be treated pcrturbativcly, 
while at large Vg, strong correlation effects are present [24]. The most studied example 
of Eq. (9) is certainly the electron liquid with pure Rashba spin-orbit coupling (n = 1). 
However, this is also a special case in which several many-body effects are left essentially 
unaffected by the presence of spin-orbit coupling. As discussed below, this behavior can 
be contrasted to hole systems with a dominant spin-orbit coupling that is quadratic or 
cubic in momentum (n = 2,3), for which special cancellations that are present for n — 1 
do not occur. 

3'1. Exchange- correlation energy. - One of the most fundamental properties of the 
Hamiltonian, Eq. (9), is the ground-state energy. By considering the many-body ground 
state with Fermi surfaces determined by x in Eq. (6), we can define the corresponding 
exchange-correlation energy £'a;c(5, J's, x) follows [45] (in effective Rydbergs), 

nn F( ^ 2g v/|i + x|"+^-v/|i-xl"+^ , ^ I . 

Here, the first two terms are, respectively, the non-interacting kinetic and spin-orbit 
energies. Note that Eq. (11) is written at fixed x and to find the ground-state energy 
it is necessary to further minimize £{g,rs,x) with respect to x- For example, Xo(3) of 
Eq. (8) is obtained by minimizing the non-interacting part of Eq. (11), i.e., that without 
fxc [39]. (1) 

As mentioned above, the n — 1 case is the most-commonly studied case. A general 
result for £[g^rs^x) in this case has been obtained from an analysis of the perturbative 
expansion of £xc{g,Ts,x)i to all orders in the Coulomb interaction parameter r^. By 
Taylor expanding a generic energy diagram in g, x, it has been found that the quadratic 
term is always oc {g — xY' [33]. Performing a similar expansion on the non-interacting 
energy of Eq. (11) gives 

(12) 8{g, r„ x) = ^ + EUts) - ^- + { \ + c] {g - xf + ■ ■ ■ . 

rj rj \ri J 

with C an unspecified constant which depends only on r^ , while the exchange-correlation 
energy without spin-orbit coupling, Sxci^s), is a well-studied quantity [105, 24, 106, 107, 
108, 109]. Provided the constant satisfies C > — I/t's (which is always true at sufficiently 
high density), Eq. (12) shows that the total energy is minimized by the same value 
of the chirality xo — 5 as for the non- interacting system, see Eq. (8). An immediate 
consequence of setting x — 5 is that the quadratic corrections to the ground-state energy 
are generally absent for n = 1. 



(^) For n — 1 and x ^ 1 one should replace 1 + with 2x in the non- interacting kinetic energy. 
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Fig. 2. - Left panel (ftom Ref. [45]): Quantum Monte Carlo energies (dots) for the electron 
liquid at Ts = 10 including the Rashba spin-orbit couplings [44]. The solid lines are obtained 
using the exchange-correlation energy in the absence of spin-orbit coupling. Right panel (from 
Ref. [39]): exchange energy at rs = 1 for eigenstates of the non-interacting Hamiltonian (4), 
parameterized by the chirality X- 



An explicit evaluation of Excirs^g^x) becomes possible in the high-density limit, in 
terms of a controlled expansion in powers of the small parameter Tg- In that regime, the 
leading contribution is given by the exchange energy [45] : 

(13) £,,(.,,5,^)^ --_ + __ = --_ + __ (\n- + - 

which is in agreement with Eq. (12), for C = 0. Higher-order correlations have also been 
examined, confirming the general form of Eq. (12) and that the coefficient C is generally 
non-zero. The second-order correlation energy yields a density-independent prefactor 
which is C = for the exchange contribution, but C = —1/4 for the direct term [45]. 
The next contribution to the correlation energy and the corresponding expression for C, 
from the infinite resummation of ring diagrams, have also been studied explicitly [45]. 
However, these higher-order terms do not modify the leading result for the ground- 
state energy, obtained from Eq. (13) by setting x — 9- The possibility of non-analytic 
corrections to the energy oc g'^ln^ g, which are still compatible with Eq. (12), has been 
discussed in Ref. [48] . 

A striking confirmation of the small effect of spin-orbit coupling extending to the non- 
perturbative regime comes from an analysis of the exchange-correlation energy obtained 
by making use of the diffusion Quantum Monte Carlo method [44]. Fig. 2 reproduces 
some of the numerical data (dots), for representative values of the parameters g,rs,Xj 
compared with the simple approximation (solid lines) of setting £xc{g, fsiX) — £xc{i's) in 
Eq. (11) [45]. As seen in Fig. 2, this approximation is in excellent agreement with the 
Quantum Monte Carlo results, even in the strongly correlated regime > 1 and at large 
values of the spin-orbit coupling parameters g^x^ 1- 
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These arguments thus indicate that the exchange-correlation energy is always essen- 
tially unchanged at finite g for n = 1. On the other hand, this behavior is peculiar to 
the paramagnetic state with n ^ 1 Rashba spin-orbit coupling. In contrast, a significant 
interplay between Coulomb interaction and Rashba spin-orbit coupling exists for broken- 
symmetry states [40, 41, 42, 43] and in the presence of a magnetic field (e.g., for the spin 
susceptibility [47, 34, 48]). Furthermore, non-lincar-in- momentum spin-orbit coupling 
terms have been studied where larger modifications of the exchange-correlation energy 
can be found. In particular, the exchange energy has been calculated in Ref. [39] and is 
plotted in Fig. 2 (at rg = 1). In contrast to the n = 1 result, the n = 2, 3 curves have a 
much stronger dependence on x- In fact, a small-^ expansion gives 



While the ^ correction vanishes for n = 1, in agreement with Eq. (13), a finite result 
is obtained with n = 2,3. 

Furthermore, we can compare the n = 2, 3 exchange energies to the usual result for a 
spin-polarized system, which is plotted in Fig. 2 as the n = curve(^). As is known, the 
exchange energy favors spin-polarized states and the n = curve is thus a decreasing 
function of spin polarization. On the other hand, the exchange energy with n = 2, 3 
has the opposite effect of suppressing the subband population difference. In a realistic 
situation, both spin-orbit coupling and the Zeeman term coexist. On this basis, it was 
proposed in Ref. [39] that a competition between the opposing effects of spin polarization 
and spin-orbit couping (induced by the magnetic field itself) could be responsible for the 
apparent suppression of interaction effects on the subband populations, observed for 
dilute hole systems in strong magnetic fields [25] . 

3'2. Spin-dependent quasiparticle velocity. - Another remarkable effect of electron 
interactions is the renormalization of quasiparticle properties. The combined effects of 
Coulomb interactions and linear- in-momentum (Rashba and/or Dresselhaus) spin-orbit 
coupling have been studied in Refs. [29, 30, 31, 32, 34, 28] with a variety of approximation 
schemes. In particular, perturbative expressions for the lifetime, effective mass, and 
spectral weight Z, have been obtained in Ref. [30] within the RPA approximation. As in 
our previous discussion of the exchange-correlation energy, it was found in Ref. [30] that 
all leading corrections (linear in g) are absent. This result is valid more generally since a 
non-perturbative argument analogous to that leading to Eq. (12) holds for quasiparticles 
as weh [33, 28]. 

In contrast to the linear-in-momentum (n — 1) case, there is no fortuitous cancellation 
when n = 2, 3, as appropriate for hole systems. This observation motivates a study of 
quasiparticles with n = 2,3 [28]. We focus here on quasiparticle propagation [28]. By 



(^) For n — 0, the parameter x has the meaning of a fractional spin polarization. 



(14) 
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(a) electrons (11= 1 ) 




t=0 t=r 



Fig. 3. - Left panel (adapted from Ref. [28]): Schematic representation of the motion of an 
initially unpolarized wavepacket for a fixed time r. (a) refers to n = 1 and (b) to n = 2, 3. 
At high density, the Coulomb interaction renormalizes the average velocity as well as the ve- 
locity difference between spin components, with respect to the non-interacting values (dashed 
wavepackets). Right panel: Schematics of the effect of repopulation on the spin-dependent ve- 
locities. The Coulomb interaction induces energy shifts of the non-interacting spin branches 
(dashed) which both reduce the chirality x a-nd contribute to the increase (decrease) of the non- 
interacting velocity v^. (ii° ) [28]. In addition to the energy shift, the slope of the non-interacting 
energy dispersion is modified as well [110, 28], but this effect is not shown here. 



considering the group velocity at the Fermi energy, the difference between linear and non- 
linear-in-momentum spin-orbit couplings becomes apparent even for a non-interacting 
system. At small rs,g, we have 

(15) !i:^l±^(„_l), 

with Vf = hkF /rn. Equation (15) implies that the motion of a spin-unpolarized wavepacket 
at the Fermi energy will be qualitatively different for the two cases {n — \ and n = 2, 3), 
as schematically illustrated in Fig. 3. If n = 1 we have v^ = v^_ and the wavepacket will 
remain spin-unpolarized for all time. In contrast, if n ~ 2,3, the two spin components 
have different velocities and will spatially separate. 

As is well-known [24], many-body effects affect the quasiparticle dispersion, modifying 
Eq. (15) for an interacting system. We thus introduce the following general notation [28] 

(16) ^ ~ 1 + 5v{r,) + 5vl{g) + Sv±{rs,g), 
Vf 

where 6v(rs) is the known spin-independent contribution from the Coulomb interaction 
at p = 0, which has been studied extensively in the literature [110, 111, 24, 108]. The 
non-interacting effect of the spin-orbit coupling is Svj_{g) ~ ±(n— 1)5/2, see Eq. (15), and 
the last term describes the nontrivial interplay between spin-orbit coupling and Coulomb 
interactions. The term Sv±{rs,g) has been studied in Ref. [28] within the context of the 
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screened Hartree-Fock approximation [110, 29] which gives, for n = 1, 



(17) Sv±{rs,g)''^''- 



16tt 



2 



The result of Eq. (17) is computed to lowest order in Ts, as an expansion in powers of 
small g. As seen in that equation, the leading term of Sv±{rs,g) contains a factor 
and is independent of the spin branch. The result is in agreement with the RPA analysis 
of Ref. [30] [only the 0{g^ In g) term was obtained there]. This agreement supports 
using the much simpler screened Hartree-Fock theory in the high-density regime. In 
Eq. (17), the ~ g^ term introduces a velocity difference between the two subbands, 
which is a qualitatively new effect introduced by the Coulomb interaction [28]. In fact, 
f+ = w*! exactly for the non-interacting system. However, this finite velocity difference 
only arises as a subleading correction. As mentioned previously, the presence of non- 
linear-in-momentum spin-orbit coupling terms (n = 2, 3) induce much more dramatic 
effects than suggested by Eq. (17). Here we cite the corresponding result for n = 3 to 
leading order in rs,g [28] : 

2\/2 

(18) (5w±(rs,5) ~ ±— — r^g, n = 3. 

OTT 

A similar expression is found for n ~ 2 and in both cases 5v± (r^ , g) ^ g and has opposite 
sign in the two chiral subbands. 

To better understand the velocity renormalization, here we consider the specific con- 
tribution to 6v±{rs,g) related to the repopulation effect mentioned after Eq. (14). Min- 
imization of the total energy, Eq. (11), with respect to x gives [39]: 

at small g and 7'^, when S^c can be approximated as in Eq. (14). While the effect of 
interactions is absent for n = 1 (when we have x — 5 as in the non- interacting case), the 
factor in parenthesis reduces the value of x for n = 2,3. The same result, Eq. (19), is 
obtained from the quasiparticle analysis of Ref. [28] , and is due to unequal energy shifts of 
the two spin branches, induced by the Coulomb interaction. As schematically illustrated 
in Fig. 3, the reduced value of k- corresponds to a decrease of the non-interacting velocity 
in the " — " (lower energy) subband, and the opposite is true for the " -I- " subband. The 
sign of this effect is in agreement with Eq. (18). However, the final result also takes 
into account the wavevector dependence of the interacting self-energy, which gives an 
additional spin-dependent contribution to 6v±{rs,g) [28]. As it turns out, both effects 
of the Coulomb interaction have the same sign and contribute to enhance the velocity 
difference. 

By referring to Fig. 3 and assuming n = 2, 3, we see how the non-interacting re- 
sult is modified at high density, within the screened Hartree-Fock approximation. For 
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a given travel time r, the total travel distance and the separation between the two spin 
components are larger than the non- interacting value (A > S). The difference in Fermi 
velocity of the two spin branches could be experimentally probed with transport mea- 
surements of the effective mass [111, 26] and Raman scattering experiments [112], as well 
as by direct optical imaging of the wave-packet separation. In fact, for a travel distance 
of order ^ /im. the travel time is 1 — 50 ps (depending on density) and typical values 
of 7 for holes give spatial separations between the two spin components of a few hun- 
dred nm (~ 200 nm for g = 0.1 and n = 3), which is within reach of Faraday rotation 
measurements [113, 114]. 

4. Spin-echo decay in single-hole quantum dots 

Spintronics research often involves studying the motion of spins in two-dimensional 
systems (the subject of the last section). Carriers that are free to move in two dimensions 
can, however, have their spins efficiently randomized through momentum scattering in 
combination with spin-orbit interactions [115]. To achieve long spin coherence times, 
it is typically necessary to quench the orbital momentum by confining carriers to zero- 
dimensional quantum dots. After freezing the orbital degrees of freedom, spin-orbit 
interactions may still limit hole-spin relaxation, but become an ineffective dephasing 
mechanism [51, 116]. Although there are still situations where dephasing can be induced 
through spin-orbit coupling for quantum-dot-confined hole spins [117], these mechanisms 
can be suppressed through stronger two-dimensional confinement. 

Even in the limit of extreme confinement, decoherence can still proceed through a 
small residual spin-spin coupling such as the hyperfine interaction between an electron 
or hole spin and surrounding nuclear spins in the lattice. For a heavy-hole spin in a quasi 
two-dimensional quantum dot with negligible strain, the hyperfine interaction iJhf with 
surrounding nuclei is approximately Ising-like [57, 118, 119, 120]: 



Setting h ~ 1, h = "^i^Aklk is the Overhauscr operator, S = cr/2 is a pseudospin- 
1/2 operator in the two-dimensional (J^ = ±3/2) heavy-hole subspace, and 1^ is the 
nuclear spin at site k. Ak is the hyperfine coupling constant for the k^^ nuclear site and 
is given by = A*'=z;o|V'(rfc)|^, with the hyperfine constant for isotope ik at site 
k, vq the volume occupied by a single nuclear spin, and ip{vk) the heavy-hole envelope 
wavefunction [57, 121]. 

We consider an in-plane applied magnetic field B = Bx (see Fig. 4), which results in 
the total Hamiltonian 



(20) 




k 



(21) 



k 



where Hz gives the hole- and nuclear- Zeeman interactions, and 7^ is the gyromagnetic 
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Fig. 4. - Quantum-dot geometry, with Overhauser field h and magnetic field B = Bx. For 
unstrained and flat (d ^ L) quantum dots, the in-plane hole gyromagnetic ratio is small, g± ~ 
[122], and the Overhauser field lies along the [001] growth direction z, hx,y — [57, 118, 119]. 



ratio of the isotope at site k having total spin For an InGaAs quantum dot, we 
take into account two isotopes of Ga (^^Ga and ^^Ga) with significantly different 7, and 
assume the same 7 for the two fn isotopes (^^■^In and ^^^In). The hole gyromagnetic 
ratio is 7h = g^fJ-B, with g± the in-plane factor for a dot with growth axis along [001] 
and fiB the Bohr magneton. 

As the hole spin evolves under the action of the Hamiltonian in Eq. (21), fluctuations 
in the Overhauser field induce fast spin dephasing. This dephasing can be reversed 
with a spin echo if the fluctuations are effectively static on the timescale of spin decoher- 
ence. We consider an echo scheme consisting of 7r-pulses about x, which would reverse 
the effects of static fluctuations in effective fields along y or z, due to, e.g., the hyperfine 
interaction in Eq. (20) (see, e.g.. Fig. 5). 

To calculate the effects of spin echo on hole-spin decay, we move into the interaction 




Fig. 5. - Echo sequence with 7r-rotations about x, which reverse dephasing from static fluctu- 
ations in an effective field Suj^ along z due to the hyperfine interaction. The magnetic field is 
B = Bx and in this figure we assume ~ for simplicity such that there is no precession of 
the hole spin about x. 
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picture with respect to Hz and define the interaction picture Hamihonian H{t) as 

(22) H{t) = e^«-*i/Me-'^-* = h.{t)S,{t), 

(23) with {t) =^Ak [II cos (7., Bt) - sin (7,, Bt)] 

k 

(24) and S^{t) = [S^ cos{-fuBt) - Sy sin(7HBt)]. 



Accounting for a 7r-pulse at time r about an axis a (a = x, y, z), followed by a second 
TT-pulse to return the spin to its original orientation, as illustrated in Fig. 5 for a ^ x, 
the echo Hamiltonian is given by 



(25) ffe(t) = 



H{t) < t < r, 

aaH{t)aa T < i < 2r. 



From Eqs. (22) - (25), we see that if a = x, then Hc{t) = -H{t) from t r to t = 2t. 
Provided H[t) is approximately static on this timescale, the 7r-pulse will reverse the 
evolution undergone by the system up until t = t and refocus the hole spin signal at 
t = 2t. Due to the anisotropy of the heavy-hole hyperfine interaction, the decay process 
is highly dependent on the chosen geometry; e.g. if a — y or a = z, the echo docs not 
reverse the sign of H{t) from t = t to t = 2r. For this reason, we will focus on the case 
with the magnetic field along x and 7r-pulses performed about x. 



4'1. Hole spin decoherence. - The hole-spin coherence dynamics can be written as 



(26) (^a(2T)) = (;7t(2r)5„(2r)t/(2T)), 

where O denotes the operator O in the interaction picture. t/(2T), the time-evolution 
operator for 7r-pulses applied at times i = t, 2t, is given by 

(27) C7(2r) = re-'-^o'''*-^'^(*\ 



with T the usual time-ordering operator. The Magnus expansion, an average-Hamiltonian 
theory typically applied to periodic and rapidly oscillating systems [123], allows us to 
rewrite the time-evolution operator as a series expansion, 



(28) 



U{2t) = e-*^-(2r) ^ g-«E^o^'"(2r)^ 
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where all orders can be found by expanding Eq. (27) in a Dyson series and applying the 
Baker-Campbell-Hausdorff formula. The first few terms are 

(29) f H{ti)dti, 
Jo 

(30) jf(i)(t)=-^ r r [H{t2),H{u)]dhdh, 

^ Jo Jo 

(3f) H(^^\t)^-l f r f\H{t3)[H{t2),H{h)]] + [H{h)[H{t2),H{t3)]]dhdt2dt3. 

D Jo Jo Jo 

Each higher-order term in the Magnus expansion contains one additional integral over 
time, such that higher-order oscillating terms are suppressed by a factor of order ||i/||/a;, 
with uj the typical oscillation frequency. In addition to bounded oscillatory terms, higher 
orders in the Magnus expansion typically also contain terms that grow in time. To guar- 
antee that these growing terms do not cause a significant modification of the dynamics 
evaluated here, we can only apply the Magnus expansion up to a timescale where these 
contributions reach ^ 1. The precise conditions for validity of the Magnus expansion in 
this case are given in Ref . [69] . 

Wc can now express the time-evolved spin components Sai^r), where a — x,y, z, as 

(32) (SMt)) = (e*^-(2r)^^(2T)) 
by defining the Liouvillian LM(i) [69] as 

(33) Lm{2t)0=[Hm{2t),0]. 

Here, the expectation value of a given operator O is defined as (O) = TrjOp}. We make 
the assumption that the initial state p = ps® pi is a tensor product of the hole-spin (ps) 
and nuclcar-spin (p/) density matrices. We further assume that the nuclear spins are 
in an infinite-temperature thermal state. In this limit, for N ^ 1 uncorrelated nuclear 
spins, the central-limit theorem allows for the Gaussian approximation: 

(34) (e'^"(2-)^«) ^ ^exp{ - i(4i(2r))^}5„^^, 

where, for any operator Os acting in the hole-spin space, 

(35) {LUt))iOs = Tri{iLl,it)Os)pi}. 

In a sufficiently large magnetic field B, rapid oscillations in H{t) will allow us to truncate 
the Magnus expansion at leading order, giving 

(36) Lm{2t)Os ^ L^°H2t)Os = [H^^H^t), Os]. 
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Fig. 6. - Exact (solid lines) and approximate [dotted lines, from Eq. (37)] solutions for various 
magnetic field values B. Parameters used for these plots were: {Sx{0)) = 1/2, N — 10*, 7h ~ 
5±MB ~ 0, and values for 7^ can be found in Table I of Ref. [121]. We have further assumed 
isotopic abundances and nuclear spin appropriate for a Ina;Gai_iAs quantum dot with indium 
doping X = 0.5, as given in Ref. [69]. 



4' 2. Vanishing g± limit. ~ In the limit 7h = 5_lMb = 0, an exact analytical solution 
can be found for {Sa{t)) [69]. This exact solution can be used to verify the validity 
of approximations applied here and will provide a good description of the spin-echo 
dynamics whenever 7h < 7;. In GaAs, this condition corresponds to g± < 10"'^ {g± < 
5 X 10""^ has been reported in 2D wells [122]). The exact evolution of (Sxi^r)) is shown in 
Fig. 6 for a range of B, revealing the transition between two distinct regimes: A low-_B 
regime resulting in complete decay, and a high-i? motional-averaging regime in which 
decay is bounded. In the low-i? regime (top two panels), the decay time decreases with 
increasing magnetic field, whereas for high B, the hole-spin coherence improves as we 
find periodic complete refocusing of the spin signal (bottom right panel). 

To gain physical intuition for this surprising behavior, we calculate the approximate 
dynamics by setting 7h = in Eq. (34) and keep only the leading term in the Magnus 
expansion, giving [69] 



AV,A^I,{I, + 1) . 4 ^1^BT 

-, 7^:; sui ' 

3iV(7»S)2 



where the index i sums over the different nuclear species, fj is the abundance of the i^"^ 
species, and A = Ak is the total hyperfine coupling constant averaged over the nuclear 
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species. Eq. (37) is plotted in Fig. 6 along with the exact solution. From this figure, we 
see that the approximate solution agrees almost perfectly with the exact solution in the 
chosen regime. 

In the \ow-B regime, where B <C A/{'-fi\/N) (top-left panel in Fig. 6), a short-time 
expansion of Eq. (37) gives 



When B increases, Eq. (39) shows that the initial decay time tq decreases, in agreement 
with the top two panels of Fig. 6. This behavior is contrary to the usual expectation 
for electron spins (for electrons, the coherence time normally increases with increasing 
B [124]). Here, a decrease in the decay time with increasing B is due to dynamical 
fluctuations in induced by rapidly precessing nuclear spins. Because the nuclear-spin 
system fluctuates in time, decay due to this dynamical environment can no longer be 
completely reversed by a simple Hahn echo scheme. 

At large magnetic fields, when B > A/{'-fi^/N), we sec from Eq. (37) and Fig. 6 (bot- 
tom right panel) that the hole spin signal no longer decays to zero. In fact, according to 
Eq. (37), the oscillation amplitude in the exponent decreases as 1/B^. In this regime, 
the nuclear spins precess rapidly enough to average out the fluctuations in hz] this mo- 
tional averaging accounts for the bounded decay as well as the suppression in oscillation 
amplitude of the echo signal. Between the two regimes, we find partial periodic recovery 
of the hole-spin signal. 

While the discussion above addresses the g± = case, we predict the same motional 
averaging effect to occur when 7h 3> 7i, at even lower applied field values [69]. The 
main qualitative distinction from the g± = case is an additional beating in the echo 
envelope function introduced by the hole Zeeman frequency. These results demonstrate 
that hole-spin dephasing due to hyperfine couplings can be completely suppressed at 
moderate magnetic fields, in contrast with the case of electron spins. 

To conclude this section, we mention a number of other interesting topics related to 
hole spin coherence and manipulation. First, the spin-echo decay observed in Ref. [62] 
is believed to be primarily caused by electric-ficld-induccd noise, instead of hyperfine 
coupling to nuclear spins. In light of this result, in Ref. [69] we have developed a simple 
phenomenological theory of spin decay brought on by fluctuations in the hole spin pre- 
cession frequency. Furthermore, new schemes for potentially improved single-hole spin 
manipulation and control have been proposed very recently as well. One such method 
suggests coherently manipulating spin states through the Dresselhaus spin-orbit interac- 
tion, where the heavy-hole spin would be transported around closed loops using static 
applied voltages [125]. Another proposal points to the fact that a hole spin's precession 
can be altered through electric-fleld- induced g-tensor modulations, resulting in faster 
predicted spin manipulation times for holes than for electrons [126]. A quantum-dot 



(39) 



(38) 



(5.(2r)) ^ (5.(0)) (1 - (r/ro)4) ^ (5.(0))e-(-/-'')^ 

1 \^l^^h^A)^ hih + l)Y^^^ 
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molecule has also been proposed as a potentially scalable qubit architecture which could 
allow for enhanced wavelength tunability and qubit rotation fidelity [127, 128, 129, 130]. 
In addition to the latest advances in maintaining and controlling the hole spin, current 
and future prospects are paving the way for establishing hole spins in quantum dots as 
a robust qubit implementation. 

5. Conclusions and open questions 

In this Chapter we have reviewed recent work on heavy holes in III-V semiconductors, 
by focusing on many-body phenomena due to electron-interactions (for two-dimensional 
systems) and the hypcrfine interactions with the nuclear bath (for single-hole quantum 
dots). In both cases, we have tried to highlight the unique behavior of heavy holes, 
which originates from the total angular momentum = ±3/2 of the effective spin. For 
two-dimensional liquids we have discussed how typical non-linear-in-momentum spin- 
orbit couplings for holes can have a stronger interplay with the Coulomb repulsion. In 
particular, the significant difference in group-velocities for holes with strong spin-orbit 
coupling implies that the spatial separation of spin components is much easier to achieve 
than for electron carriers. This velocity difference could be observed with optical tech- 
niques including Raman scattering and direct time-resolved Faraday-rotation detection 
of hole-wavepacket motion, in addition to transport measurements of the effective mass. 

Other properties of the hole liquid will likely show similarly interesting behavior. 
Future studies might include a complete characterization of the quasiparticle param- 
eters (lifetime, spectral weight) and compressibility. The challenging regime of dilute 
hole liquids with large interaction parameters ~ 6 — 12 calls for further progress in 
non-perturbative approaches, e.g., a suitable extensions of the Jastrow- factor variational 
wavefunctions and the Quantum Monte Carlo method. 

On the other hand, single holes in quantum dots represent an emerging platform 
for the implementation of spin qubits and are attracting a growing interest, as is clear 
from the fiourishing experimental progress of the field. Following the first realization 
of the Hahn spin echo in this system [62], we have described a theoretical analysis of 
spin-echo decay based on exact results and a systematic method of solution in terms of a 
Magnus expansion [69] . We expect this approach will be a useful tool for further studies 
of spin decoherence, in this or other related systems. While electric noise seems to limit 
the spin-echo decay in current experiments [62] , the prediction of a motional-averaging 
regime at high magnetic field will be especially interesting for a future generation of 
devices, dominated by hyperfine decoherence. The presence of such a regime strengthens 
the case that holes will be able to realize long-lived spin qubits, amenable to fast optical 
manipulation. 

* * * 
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